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Abstract 



We construct the scaling limit of the easy axis XXZ chain. This limit is a 
subtle combination of approaching the isotropic point, and letting the lattice 
spacing to zero to obtain a continuous model with a finite mass gap. We give 
the energy difference between the two lowest energy states (the two 'vacua') 
and analyze the structure of the excitation spectrum of the limiting model. 
We find, that the excitations form two sets corresponding to the two vacua. In 
both sets the dressed particles are described by Bethe Ansatz like equations 
(higher level Bethe Ansatz), and the two sets can be distinguished through 
a parameter entering into these secular equations. The degenerations in the 
spectrum can be interpreted as originating from an SU(2) symmetry of the 
dressed particles. The two particle scattering matrices obtained from the 
secular equations are consistent with this symmetry, and they differ in an 
overall sign in the two sectors. 

PACS numbers: 75.10.D, 05.50.+q, ll.10.Kk, 11.15.Tk 
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I. INTRODUCTION 



The XXZ Heisenberg chain is defined by the Hamiltonian 

V N 

Hxxz = - E + S y n S y n+1 ) + p (s z n S z n+1 - 1/4)) . (1.1) 

a n=l 



The Hilbert space is the tensor product of N spaces furnishing the doublet representation 
of SU(2). S%' y ' z are the spin operators acting on the nth site and in (|1.1| ) the (iV+l)th site 
is identified with the first one. We consider the easy axis region defined by 

p = cosh 7 > 1 . (1.2) 



The V/a factor with a being the lattice spacing and V = 2 /n is included for the sake of 
proper normalization. 

The Hamiltonian (|1.1|) can be diagonalized by Bethe Ansatz (BA) |l|-f||, the BA equations 
(BAE) were analyzed in [§-0]. Due to these studies it is well known by now, that the 
chain with even number of sites has two states in which there are no free parameters ||. 
These are the ground state and the first excited state, and as the difference in their energy 
disappears exponentially fast with the length of the chain N 0,0], they are often referred 
to as the two ground states, although it is more precise to call them the two vacua. The 
other excited states are described in terms of a few parameters characterizing some kind of 
dressed particles. These posses a gap, and their parameters satisfy BA type equations, the 
so called higher level Bethe Ansatz equations (HLBAE) [^>J. The excited states form two 
groups which can be associated with the two ground states |J. 

The model has two relativistic continuum limits ||. One is constructed by setting 7 = 
first and taking the a —>■ (together with iV — > 00 but Na = L) continuum limit afterwards 
resulting a massless SU(2) conformal field theory. To obtain the other relativistic limit 
one has to perform the 7 — > and a — > simultaneously keeping the gap in the properly 
normalized spectrum constant. The aim of the present work is to study the details of this 
second limit, and analyze the limiting theory. Our results are as follows. 



The limiting theory is a massive relativistic theory with a mass M, if the continuum 
limit is performed in such a way, that 



7T 2 



7^0, a^O, -expj-— |^M. (1.3) 



As the length of the chain L = Na is kept finite, also the number of sites must be 
adjusted: 
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LM (ir 2 ) 

iV=— expj-j^oc. (1.4) 

Our expression for the mass M differs from that found in in a prefactor of the 
exponential. 

In the above limit the difference in the energies of the two 'ground states' (physical 
vacua) stays finite but exponentially small in the length of the system: 



A* = ^-. (1-5) 

The excited states of the system can be described in terms of excitations (dressed 
particles). Each particle is characterized by a rapidity we denote by The energy 
and momentum contribution of the particles are the sums of the contributions of the 
individual particles 

E-E = Y,£($), P = EpW' (1-6) 

with 

e ($) = Mcosh^, p{d) = M sinh <& . (1.7) 
The rapidities of the particles must satisfy a set of BA type equations: 

Lp(V h ) = 2n(I h + J ) - £ J + £ 2 tan- 1 *L«i , 

i \ 71 J n / 2 

X h = - K \ 1 ' (odl) , J = or 1/2, (1.8a) 



2 V 2 



and 



"W «■ _,9. n ( K ) K _ Ka 

V 2tan- 1 ^-f^ = 27TX + V 2tan~ 1 , 
h ir/2 p it 

( n(n) - n(<&) + l\ , , . . . 

J«=[ o (° dl )- (l-8b) 



Here n(i?) is the number of particles, the set of variables k is needed to describe the 
internal symmetry of the states, their number n(n) obeys ri(i?) — 2n{K) > 0, and 

m = -in , . i.9) 



i 



r(| + ii)r(i-zf) ' 
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• The excited states form SU(2) multiplets in which the energy is the same. Each mul- 
tiplet is characterized by one solution of the equations (|1.8Q . The number of states be- 
longing to one multiplet is — 2n(«)) + 1. Within a multiplet the states are labeled 
by the value of the S z taking the values n('d)/2—n(K), n{d)/2— n{n)— 1, — n{d)/2+ 
7i(k). This indicates, that the particles have SU(2) symmetry or a symmetry produc- 
ing the same multiplet structure, and the z component of the spin connected to this 
symmetry coincides with the z component of the real spins building up the original 
chain. 

• The solutions of ( |1.8|) form two groups, one for I = and one for I = 1/2. We 
argue, that the two sets of excited states are the excitations of the two vacua. In both 
sets the S z = states, including the vacua, are eigenstates of flipping all spins, but 
with different eigenvalues. In the case of the two particle excitations the symmetry of 
the singlets is the same as that of the corresponding vacuum, and that of the triplets 
is the opposite. This is the same structure as found in the XXX chain, where it is 
a consequence of the SU(2) symmetry of the Hamiltonian (and this structure can be 
interpreted in terms of SU(2) with a modified coproduct, or equivalently in terms of a 
g-deformed SU(2) at q — —1). 

• The two particle scattering matrix can be given up to an overall phase as 

S { At» = - exp {i (ar/„ + ( v)) } ( P « + ^S A ) ' (L10) 

where P tr and P s are the projectors on the triplet resp. singlet subspace of the two 
spins. 



All these shows, that the scaling limit (SL) of the XXZ yields a massive relativistic theory 
(in duplicate) which is of the same structure as the massive sector of the theory obtained 
through the SL of the attractive Hubbard chain and identified as a regularization of the 
SU(2) symmetric chiral invariant Gross-Neveu (CGN) model [|K],|ll|]. It is widely accepted, 
that the XYZ chain is a lattice regularization of the sine-Gordon (SG) theory or the massive 
Thirring model (MTM) |12| (as these latter two are equivalent |13|). It is also known, that 



the MTM/SG theory at a special value of the coupling is equivalent to the CGN model 
(up to a free massless boson field), more over this value of the coupling corresponds to 
the antiferromagnetic (our p=l) point of the XYZ chain. This way recovering the massive 
sector of the CGN model is not unexpected. It is remarkable, however, that the SL of 
the XXZ chain is taken through couplings not corresponding to stable MTM/SG theories. 
Our detailed analysis also rises some questions connected with the mass-formula and the 
symmetry of the vacua. 

The ( lOp mass-formula is not of the form found in the case of the Hubbard chain (a pref- 
actor ^7 is 'missing'), and it defines a /5-function different from the one expected based on 
perturbation calculations. This may be connected with the fact, that the p=l point in the 
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space of the couplings of the XYZ model is singular (in the sense, that certain quantities like 
the ground-state energy of the XXZ chain are singular in this point |§), thus approaching 
this point from different directions may lead to different results. In our SL the isotropic 



point is approached along the XXZ line while the limit taken in |L2| involves the XY-type 
anisotropy. 

Another interesting point is the existence of two very similar but not completely identical 
limiting theories. The two theories are distinguished by the symmetry of the vacuum with 
respect to reversing all spins. One has to note, however, that measuring this symmetry in 
the continuum limit may encounter difficulties, as the SL of the corresponding operator can 
not be constructed directly. This symmetry has not been studied in other representations of 
the CGN model although it may be present also in other cases. In the SL of the half-filled 
Hubbard chain it is even stronger in the sense that the vacuum is a singlet of both the spins 
and the isospins, and it is symmetric under both the spin and isospin reversal. 

The paper is organized as follows. In Section [H| we summarize those properties of the XXZ 
chain we need to construct its scaling limit and we also propose a formula to calculate the 
symmetry with respect to flipping all spins. In Section |TXT| the SL is constructed and the 
properties of the limiting theory are analyzed. Specially we construct the limiting process 
resulting in a relativistic dispersion (|III A| ), analyze the vacua in this limit (|IIIB|) , derive 
the secular equations of the limiting theory ( [111 (J] ), we give the SU(2) multiplet structure 
of the eigenstates (|IIID| ), derive the two particle 5*- matrix (|IIIE| ), we study the reflection 
symmetry of the two particle S z = states (|IIIF1 ), and we also argue, that the boundary 
condition obeyed by the particles of the limiting theory is connected to the parity properties 
of N kept unchanged in the N — > oo limit ( [III Gj) . The more technical details are collected 
in appendices. The properties of the elliptic functions we use are listed in Appendix [A]. The 
behavior of the two vacua are calculated in Appendix [B[ and in Appendix we check if 
the the approximations used to derive the HLBAE of the lattice model remain valid also 
in the SL. In Appendix [D] we describe the structure of the solutions of the HLBAE and 
give the classes of solutions which become degenerate in the SL. Finally some details of our 
numerical calculations are given in Appendix |E[ 



II. BETHE ANSATZ SOLUTION OF THE XXZ CHAIN 



A. The BA equations 



This section is intended to summarize the BA solution of the XXZ chain with emphasis 
on those points which are relevant to the construction of the SL. 

First let us recall some properties of the Hamiltonian ( |1 . 1|) . In the literature two versions 
of the XXZ model are commonly found which are distinguished by the relative sign of the 
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(S*S* +1 + S^Sn + i) and (S*S* +1 ) terms. The two theories are equivalent if N is even, the 
corresponding transformation on the Hilbert spaces is generated by ]J k cr| fc which relates the 
spin operators as: {S% k , S% k , S^} — » {—S^, — Sffc; As this is a symmetry of the Hamil- 
tonian, the energy spectrum does not change, but it does have an effect on the eigenvalues 
of the momentum operator which is relevant to the connection to the continuum limit. The 
momentum of every spin wave is shifted by 7r which affects all the states with an odd number 
of spin waves, including the ground state whenever N/2 is an odd integer. There is another 
useful operator, which commutes with the Hamiltonian, but whose eigenvalues are affected: 



N 



n=l 



which represents a reflection on the x-axis. (In a basis given by the products of the S z 
eigenstates of the individual spins E simply flips all the spins (the up ones down and the down 
ones up)). The symmetry connected with this operation we call spin reversal or reflection 
symmetry.) The S z = BA eigenstates of (|1.1|) are expected to be also eigenstates of Q2.1|) . 
For the sake of definiteness we have chosen the positive sign in (|1 . 1\) and note that when 
N is an integer multiple of four, the results should agree with the ones obtained from the 
opposite convention. 



The Hamiltonian (|1 . 1|) and S z can be simultaneously diagonalized by the Bethe Ansatz (BA) 
and due to the symmetry corresponding to inverting the spins it is sufficient to consider 
only S z > states. The eigenvectors are explicitly given in terms of a wave number set 
|{fci . . . fc r }), for the states with S z = N/2 - r, (r < N/2) 

| fa . . . M> = E ( E ex P \i E k v a n a + - £ 9 VaV , } I] cr~ a \F) ), ( 2.2 ) 

m<...<n r y -p [ a £ a <l3 

where the summation is over all permutations V of {1, . . . ,r} and the wave numbers k a 
together with the phase shifts \l/ Q/ g satisfy the equations 



cot ^ = p co^-cot^ 

p ;i- P )cot^cot^-(i + p) ' 



exp{iVA; a } = exp <( E • (2.4) 



The energy and the momentum of the state are 



E =-J2(cosk a -p); Q = E fc - (2-5) 
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The equation ( |2.3| ) can be solved for the phase shifts by parameterizing k^s in terms of new 
variables, v a [f| which are commonly called rapidities: 

k a = $(v a , 7 /2) , V a> i3 = ®(v a -v p ,i) , (2.6) 



with 



1 sin(2: + iS) 
$(z,6) = - In — -( 



(2.7) 



When expressing Q2.4| ) in terms of the rapidities it turns out to be useful to take the logarithm 
of the equation. In order to do this one has to define the branch cuts of $(z, 5). We choose 
5) to be a continuous function of z with $(0, 5) = tt in the strip \lmz\ < 5, while in the 
regions Imz > 5 and Imz < —5 we choose those levels in which $ — > ^f2i5 if Imz — > ±oo, 
respectively. This results in cuts running along the lines Imz = ±5, Rez < and Re^ > 7T, 
and corresponds to choosing fi — f 2 — in 0. Taking the logarithm of ( |2.4|) it becomes: 

N$(v a , 7 /2) = 2nl a + ^ $(t; Q - Vf) , 7) , (2.8) 

/3 



where the I a are half-odd integers and turn out to be very useful quantum numbers for 
characterizing the states with macroscopic number of spin waves 0]. The energy and the 
momentum in terms of the rapidities are 

e = E 7/2) , Q = E 7/2) . (2.9) 



For technical reasons it is useful to restrict the real part of the rapidities to an interval of 
length tt. We choose < Rev a < it, but note, that for our results it is important only, 
that the interval contains the point ir/2 in its interior. Note also, that the wave function, 
the energy and the momentum are independent of the particular choice of this interval and 
branch cuts of $(2, 6) (the latter modulo 27r). 



B. The ground states 



In the ground state(s) of the XXZ model with iV even, there are N/2 rapidities (S z = 0), all 
of them are real and all the Ij quantum numbers are consecutive half-odd integers satisfying 
Ij + i — Ij = — 1 for > \j (Here and in the following we denote the real rapidities 
by Xj.) For large N the distribution of A's can be well approximated by a smooth density 
<7o(A) satisfying the linear integral equation 
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A+7T 

-$'(A, 7 /2) = <x (A)- — J $'(A-A', 7 )<7 (A')<iA', (2.10) 

A 

N/2 

with iV$(A, 7 /2) = 27r(/ 1 + l/2) + ^$(A-A i ,7). (2.11) 

j 



Eg .( p. 10 ) can be solved by Fourier transformation leading to 



K (2K \ K' 7 

a (\) = —dn(—\,k) with — = -, (2.12) 

TT Z \ TT J K TT 

where dn(w) is the Jacobian elliptic function and K is the complete elliptic integral of the 
first kind with modulus k. (See Appendix |A] for notations and some properties of the elliptic 
functions.) 

By means of this density the sum over j in Q2.11| ) can be calculated, and the rapidity set can 
be reconstructed. One finds, that there are two nonequivalent sets of N/2 real A/s satisfying 
(PD with - Ij = — 1 f§. These are: 

TT „ /2TT 



X >~2K F \N 



for Af/ 2 =even and 

TT „ (2lT 



(j-Jo),^, Ij = 2I Q — 1/2 — j , j = 1,2,... N/2 (2.13a) 



Aj 2iT 



F(-^(j-l/2 + / )),fc) , Ij — 1/2 — 2/ — j , ] = 1,2,... AT/2 (2.13b) 



for A^/2 =odd, with F(w,k) being the elliptic integral of first kind. The two sets are 
distinguished by Jo taking the values of 1/2 and 0, respectively. These two states of lowest 
energy are almost degenerate. The true ground state corresponds to J = 1/2 and the 
difference between their energy is exponentially small in the length of the chain ||[7| . The 
leading term is: 

V /rP k N ^ 2 

AE = E(I = 0) - E(I = l/2) = ^^sh T i^ , (2.14) 



where k' = yl — W and \fk[ = (1 — k')/k. As A£"o is exponentially small in the length of 
the chain and - as we shall see - all the other states are separated from the lowest two by 
a finite gap, we call these the two ground states or the physical vacua. The two vacua can 
be distinguished for instance by their momentum, Q and the eigenvalue, S of the operator 

Q(I ,N) = 2vr(/ - 1/2 + N/4) (mod27r) (2.15) 
£(/o,AT) = (_i)(2/o-i+iV/2) _ (216 ) 
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C. The excited states 



The states other than the two vacua are called excited states and the ones of not too high 
energy can be described in terms of some dressed particles. These are created by leaving 
holes in the sequence of the A/s, and by introducing complex rapidities. The smooth density 
cx(A) of the real rapidities is still a good approximation but ( |2.10 ) receives contributions from 



both the holes and the complex rapidities. When one solves for cr(A) we are left with a few 
equations relating the parameters of the excitations only f|||: 

l=exp*jiv(am(^,^-0- (2.17) 

(n(0) n( X ) \ 

1 = exp % IJ2 HXa - Oh, 7/2) - E ®(X« - X/9, 7) " * > (2.18) 



(h=l 13=1 



where 



F{x, 7) = x + — sin(2mx) , (2-19) 

~^ m cosh (7m) 

and Io is again either 1/2 or 0. The positions of the holes in the real-rapidity distribution 
are denoted by 6h and their number is denoted by n{6), whose parity is the same as that of 
N. The variables \ represent the set of complex rapidities z in the following way. If for a \ 
|Im%| < 7/2, it represents a close pair 

z ± = X ±il/2 + 0(e~ Nn{z+ '~ /) ) (2.20a) 

with £l(z, 7) being a positive valued function of the order of unity (as long as 7 > 0). If 
|Imx| > 7/2, the x represents a wide root: 

z = x + isgn(Imx)7/2. (2.20b) 

As the x s are either real, or form complex conjugated pairs, the typical complex rapidity 
configurations are the 2-strings (the two z represented by a real x)i the quartets (the four 
complex rapidities represented by a complex conjugated pair with |Imx| < 7/2), and wide 
pairs (corresponding to a complex conjugated pair of x s with |Imx| > 7/2). The total 
number of the complex parameters Xa is denoted by n(x)- 

The energy and the momentum of the excited states can be expressed in terms of the 
parameters of the holes only and are given as: 

VK . . , /2K f 

a 7T V 7T 

IL 

Q = Q(I 0j N) + Y,q(0h) with g(#) = am (— 9,k)-^, (2.22) 



E = E + Y j e(6 h ) with e(0) = sinh 7 dn 9, k , (2.21) 

(¥**)-? 



h 
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with E being the ground state energy, and Q(Iq, N) is given in ( |2.15|) . The appearance of 



Q(Io, N) in the momentum suggests, that the 1$ = 1/2 resp. 1$ = excited states should be 
considered as excitations above the Iq = 1/2 resp. Iq = vacua. The relation between the 
functions e{9) and q(8) give the dispersion relation of the particles: 

V K 



e(q) = sinli7y 1 — k 2 cos 2 q . (2.23) 

The spin of the state characterized by a solution is 

S z = l -n{e)-n{ X ). (2.24) 

The states with S z < are obtained by flipping all spins, for example using the operator 
S of (|2.1| ). The S z = BA eigenstates are expected to be eigenstates of this operation: 
otherwise certain points of the spectrum were twofold degenerated (with the two states 
connected by £). This can happen accidentally but is not forced by the symmetry, as S 
has one- dimensional representations only. Recently Doikou and Nepomechi has proposed a 
formula for the parity of the S z = in the planar (p < 1) region [ 16f . Now we conjecture, 



that in the easy axis region we study, the spin reversal symmetry of the S z = states in 
analogy with their formula is given as 

N 2 N/2 

S=(-lf, with (1 = _ + (2.25) 

2 7T a=X 

We have not proved this formula, but performed numerical calculations to see this symmetry 
(Appendix [E]): solving the BA equations numerically for a number of S z = states has 
shown, that the wave function is either symmetric or antisymmetric under £, and the value 
of £ is correctly given by ( [2.25] ). In the large N limit ( 2.25 ) can be transformed into a 
simpler form using the density of the real roots f|: 

TV 2 f n ^ 1 n ^ 1 
fi = (2J + 1) + - + - E Xa ~ g E e h \ (mod2). (2.26) 



In summary, the structure of the states can be described as follows: the system has two vacua 
found in the N even system, one of them is symmetric, the other is anti-symmetric under 
( |2.1D , one of them has no (lattice) momentum, in the other this quantity is ti. Above each 
vacuum there is a set of excited states which can be considered as scattering states of dressed 
particles. The momenta of the particles are quantized through a set of equations of the BA 
type ( |2.17| - pTT8D (HLBAE). The excitation energy is the sum of the energy contributions of 
the individual particles and the momentum is the sum of the contribution of the individual 
particles added to the momentum of the corresponding vacuum. The S z = states are 
eigenstates of the operation flipping all spins, and the eigenvalue of this operation is related 
to the chain length, quantum number Iq and the rapidities describing the excitations. 
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III. THE SCALING LIMIT 



A. The relativistic spectrum 



Before constructing the scaling limit we note, that certain quantities (the momentum Q and 
the reflection E) depend on the parity of iV (and N/2) but not on its magnitude. To make 
the procedure uniquely defined we carry out the iV — > oo limit through Ns being integer 
multiples of four, in which case both in Q and £ the N can be replaced by zero, and at the 
end of the section we discuss the consequences of other possible choices. 



First, we show that ( |1.3|) defines the limit in which the dispersion relation ( |2.23| ) describes 



relativistic particles. Note that the continuum momentum (P, p) is obtained from the lattice 
momentum (Q, q) after division by the lattice spacing a: 



e 2 



^ smh7 ) 2 (l _p cos2M)K (i^) 2 |(^ )+p2 ,. (3 ,) 




Here we expanded cos(ap) about 1 and now observe that in the continuum limit (a — > 0) a 
relativistic dispersion relation emerges, provided 7 — > together with 

fVKksinh^\ , fk'\ „ , n . 

and (-J=M. (3.2) 

These requirements imply V = 2/tt and ( |1.3|) (see Appendix [A]). 

It will prove to be useful to introduce the scaled rapidity variable, •& which also makes the 
limiting procedure more transparent. It is defined in terms of the 9 and 7 as (see also (|2.12|) ): 

20 

K + ^ = K— } (3.3) 

while the energy = e(0(i?)) and momentum = -q{6{'&)) can be expressed as 

2Ksmh^k' 1 _ sm^api^d)) k' sn(i?) (3 4) 

7r 2 a dn($) a a dn(^) 

In the scaling limit d/K approaches as k — > 1 which yields 

* (3.5) 

7 

-> - cosh(^) ; p(i?) -»■ -sinh(^), (3.6) 

a a 
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in accord with k'/a — > M = constant. 

The total momentum also contains a macroscopic quantity ( j2.22j ) and can be written in the 
continuum limit as, 

p= QWL + £Msinh(# fc ). (3.7) 



In view of ( |3.7| ) the two sectors corresponding to Jo = 1/2 and Io = become infinitely sepa- 
rated in momentum space lending a strong support to interpreting the Jo = 1/2 resp. Iq = 
excited states as excitations above the corresponding (J = 1/2 resp. Iq = 0) vacua. Tech- 
nically while ( p.7|) gives the true momentum for the Iq = 1/2 sector in the a — > limit, the 
continuum momentum of the Iq = sector needs to be redefined as 

P = g(/o) + - + M sinh(^) . (3.8) 



As an alternative, we might keep both sets by redefining the lattice so that two sites form 
one elementary cell. In this case the lattice momenta and n are equivalent, and the Q(Iq) 
term can be dropped. In any case the energy and momentum measured relative to those of 
the corresponding vacuum are 

£-£ = E Mcosh (^), P = ^ Msinh (^)- ( 3 - 9 ) 

h h 



B. The vacua in the scaling limit 



The two lowest lying states are characterized purely by the density of rapidities, o"o(A) whose 
behavior in the scaling limit is crucial for our arguments. The reason is that our considera- 
tions about the excited states are based on the existence of a non-vanishing background of 
rapidities (see also Appendix 0). One can show that although cxo(A) vanishes around n/2 
(and diverges at and 7r), the density of the $'s is non- vanishing as required. The density 
£?o($) °f the rescaled rapidities (|3.3j ) can be expressed as 

Na (Jg(# + #)) d\ = L Qo {d)d$ , (3.10) 

provided dd = 2KdX/^ (here the l.h.s. is the number of A's in the interval {A, A + d\} at 
A = 7r/2 + tt$/2K, and the r.h.s. gives the number of i?'s in the interval •& + d$}), which 
leads to 



Ootf) = -LMcosh^) 



(3.11) 
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The ( PD difference in the energies of the two vacua is obtained directly by evaluating the 



scaling limit of ( p.!4| ). We have to note, however, that (|2.14D is the leading term in a large 



N expansion, and as the expansion coefficients are functions of 7 one has to see, if these 
neglected terms behave properly. In Appendix |B| we show that (|1.5| ) is the leading term of 
a large L expansion of the same quantity and the neglected terms decay faster. 

Finally we note, that the reflection symmetry is not effected by the limiting process, i.e. £ = 
±1 in the two states. 



C. Secular equations in the scaling limit 



We define the secular equations of the limiting theory as the SL of the HLBAE of the lattice 
model. This limit exists but if it is meaningful is a delicate question: deriving the HLBAE 
two approximations has been used, integrating instead of summing over the rapidities, and 
neglecting the exponentially small corrections to the close pairs. These have to be justified 
even in the SL. Replacing the sums with integrals require that the density of rapidities 
does not vanish in the scaling limit (see the previous section) more over it can be shown 
that the error introduced remains negligible. Also the corrections to the close pairs remain 
exponentially small even in the SL as it is shown in Appendix y. 

To derive the scaling limit of the HLBA we first observe, that the energy and momentum 
diverge unless all 9 scale to tt/2 as 7 — > (i.e. all $ are finite). In this case we expect also 
some of the to scale to it/2. For convenience we introduce their rescaled versions as 

K = IX {6.12) 

7 

It is straightforward to determine how the functions am , T and $ behave in the scaling 
limit: 

- (aM + K) - -) ^pftf) = Msinhtf, (3.13) 
a V 2/ 

(3.14) 




(3.15) 

\ VI / III I 

<S>(a X , — ) ^0 if A X ^finite. (3.16) 

V mj 

After substituting the above expressions to the HLBA equations, their scaling limit is ob- 
tained. The complex rapidities which remain at finite distance from tt/2 drop out from the 
equations, although they do affect S z of the state as we will see in the next section. The 
equations for the holes and the complex parameters approaching tt/2 become 
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Lp(^) = 2n(l h + J ) - £ + £ 2 tan" 1 , (3.17) 



7 a T/2 



1 /n(tf) - 2n(«)\ / t \ , 
Ih=2 [ 2 ) ( ° d 1} ' (3 - 18) 

y 2 tan" 1 l^ll = 27r J a + V 2 tan' 1 I^Jll , (3. 19) 



It is to be noted, that the quantum number I directly appears in the secular equations 
of the excited states. As we have discussed before, we interpret the two sets of solutions 
distinguished by the value of I as the excited states of the corresponding vacua. 



D. Multiplet structure of the states 



A few simple solution of the above secular equations are given and the structure of solutions 
is discussed in Appendix [D| Here we discuss the degenerations in the spectrum developing 
in the SL. As both the energy and the momentum depend on the rapidities of the holes only, 
we expect degenerate multiplets corresponding to solutions differing in their \ content only. 
As ( |3.17| ) and ( |3.19| ) contain both the i? and k sets the solutions becoming degenerate can 



differ in the x's not scaling to tt/2 (and which disappeared from the equations). 

In Appendix [D] we discuss the types of different solutions of the HLBAE fl2.18|) and relate 



them to the solutions of ( 3.19 ). We argue, that at a given number of holes, n($), each having 



the form 9h = 71^/^ + 7r/2, to any number < n(/c) < n(i?)/2 (|2.18 ) have n(i?)/2— n(/c)+l 



different solutions in which n(/c) of the x s scale to 7r/2 like \ a = 7k q /7t + n/2 with the 
common K a set being one solution of ( |3.19| ). These solutions differ in the number of x s n °t 
scaling to tt/2. It is clear than, that the number of different solutions of ( |2.17|42~T8"1) becoming 
degenerated in the SL and given by one solution of the system ( |3.17| - |3.19"D is n(i?)/2— ti(k)+1. 
The number of xs n ot scaling to tt/2 can be 0, 1, . . . , n(i?)/2— ti(k), and the corresponding 
states have spins S z = n($)/2—n(K), . . . , 1, respectively. Completing this set of states with 
the states S z = —1, —2, . . . , — (n(i?)/2 — n(«)) which are obtained by flipping all spins we 
see, that the degeneration of the states given by one solution of the equations ( p.l?H3.19| ) is 



n(i?) — 2n(tz) + 1 This corresponds exactly to SU(2) multiplets: the states are characterized 
by two quantum number I = n(i?)/2 — ti{k) and I > m > —I, and the m does not influence 
the energy. This way we may consider the particles as 'spin' 1/2 particles forming SU(2) 
eigenstates with quantumnubers I and m. (This picture will be refined when discussing the 
spin reversal symmetry.) Apparently the equations ( |3.17| ) and ( |3.19| ) give the / = m states 



directly (as usual in SU(2) BA systems), and to have the complete description of the m < I 
states one should find the x s n °t scaling to tt/2, nevertheless these parameters do not enter 
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into the spectrum. This interpretation is consistent with the expectation for the number of 
different solutions: the number of nonequivalent solutions of (|3.19| ) given by flD22|) is just 
the number of different I = n($)/2 — n(K), I = m states of ni^) spin 1/2 particles. One has 
to note, however, that while the quantumnumber m coincides with the S z of the system, 
it can not be seen, if S 2 = 1(1 + 1) or the quantity given by /(/ + 1) is different from the 
square of the total spin. Note also, that this SU(2) symmetry develops in the SL, but it is 
not present in the initial model. 



E. The two particle scattering matrix 



The two particle scattering phase shifts of the physical particles can be reconstructed from 
the secular equations. The method for this is based on the idea, that the deviations of the 
particle momenta from the free values can be interpreted as the phase shifts of the particles 
scattering on each other fll7-19[j. Consider a two particle scattering-state on a ring. If the 



momenta are p\ and P2 and the particles obey twisted boundary conditions with a twisting 
angle ip (<p = for periodic and it is 7r for the antiperiodic boundary condition), then the 
boundary condition requires Lpi + 5 12 = 2imi + (p and Lp 2 — 8x2 = 27m 2 + tp with n\ and n 2 
integers and byi being the phase shift |20|,21|. Writing the equations ( |3.17| ) in this form the 



phase shifts can be found. A weakness of this procedure is that in most of the interesting 
cases ip is not known (as the boundary conditions for the dressed particles can differ from 
those of the bare ones). The twisting could be seen directly in the equation describing 
one single particle. Since, however, the parity of the particle number is usually determined 
by that of the number of elements in the chain, the sectors with even and odd number of 
particles are disjoint in the sense, that a given system represents either one or the other, and 
one can not see, if the twisting is the same in both sectors. (A trivial example for particle 
number dependent boundary conditions is given by a spin chain where the spin waves are 
described by Fermions through a Jordan- Wigner transformation.) Now we suppose periodic 
boundary conditions, i.e. <p = 0, but we keep in mind, that our results hold up to a rapidity 
independent constant (see also 




A triplet state of two particles is described by two i?s and no ks. For such a state Eq. (|3.17|) 
yields 



^2 = +7T-27T/ , A# = #i-0 2 , (3.21) 



where the 7r comes from the parity prescription for the parameters Z^. A singlet state is 
characterized in addition to $1 and $2 by a k for which Eq. ( p\19| ) yields k = (i?x + $2)/2- 
For this case the above reasoning leads to 

/AVA /A$\ 
8 S 12 = <f) I J - 2 tan" 1 f J - 27r/ . (3.22) 
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From the phase shifts (|3.21|) and ( |3.22| ) the two particle S'-matrix can be given as 




S(A») = -exp{ i [2nI + ^[ — ))}[P tr + ^^A) , (3.23) 



where P tr and P s are the projectors on the triplet resp. singlet subspace of the two spins. 
This S'-matrix is consistent with the SU(2) symmetry of the limiting model detected in the 
degenerations. 

The scattering matrices ( |3.23|) of the two sets of solutions differ in an overall sign due to the 
appearance of Jo- To decide, if this difference is significant in the sense that it has physical 
consequences regarding the nature of the particles in the two sets of states, needs further 
considerations. 



F. Spin reversal symmetry 



Finally we have to discuss the question of the symmetry of the excited states with respect 
to £ of ( |2.1| ). The two vacua are eigenstates with eigenvalues ±1. As we have mentioned it 
already, checking numerically a number of S z = excited states of the finite lattice model 
(with finite 7) we have seen, that these states are also eigenstates of E with the eigenvalue E 
given by ( p.25| - [2~2"6| ), but on that level we could not see an easy rule for the eigenvalues. We 



have found some indications, however, that there is a correlation between the spin-structure 
of the dressed particles and the spin reversal symmetry. We have examined a number of 
solutions of different types with two holes and one close pair of the N = 20 lattice. It has 
been found that in accordance with Appendix [D| the solutions form two classes: in class 1 
the real part of the close pair is in the vicinity of (61 + #2)/ / 2 while in class 2 it is in the 
neighborhood of (#1 + 9 2 ± tt)/2. In case of Iq = 1/2 the solutions in class 1 are symmetric 
and the others, in class 2, are antisymmetric under E while just the opposite is true for 
I = (Appendix p. 



In the scaling limit the solutions of class 1 resp. class 2 satisfy (|D5|) resp. (|L)6|) representing 



singlet and triplet states. We believe that while increasing the value of N the symmetry 
properties do not change, and we conclude, that in this type of states the singlets have the 
symmetry of the vacuum, and the symmetry of the triplets is the opposite. Accepting the 
validity of ( ^25|) (and so that of flUD) we can have E also for the n{9) = 4, n(x) = 2 



states. Taking the results of Appendix [D| we find that 

E = E(/ , N, I) = (_i)(2/o-i)+iv/2-j = E (j 0j N)(-l) 1 (3.24) 



with I = 0, 1, 2, corresponding to singlet, triplet and quintuplet configurations, respectively. 
We think, this is an indication, that (|3.24|) is generally true. 
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If we perceive the reflection symmetry as the product of the symmetry of the vacuum and 
that of the spin structure of the excitations, the latter should be 

E = (-1)' (3.25) 



i.e. we have to consider now the singlets to be symmetric and the triplets anti-symmetric. 
This is just the opposite as in the case of two ordinary SU(2) spins, and resembles a q- 
deformed SU(2) structure at q — — 1 P3"| . It can be shown in general | ]241| , that defining the 
spin operators for a system of particles as 

^ = E(-1)°'~V °£ = £°i (3-26) 

(with <jj being the Pauli matrices acting on the spin labeled by j) results an SU(2) structure 
(equivalent to a g-deformed SU(2) at q = —1) in which the S z = members of the multiplets 
have the eigenvalue (|3.25|) . Particularly in the case of two spins the erf = member of the 
triplet is ^(| T!) - | IT)), while the singlet is ^(| U> + | IT)). 

We have to note, that the symmetry properties of the limiting theory are in strong analogy 
with that of an ordinary XXX Heisenberg chain. The eigenstates of this model are S 2 and 
S z eigenstates with eigenvalues 1(1 + 1) and m, respectively, with I > \m\ integers if N is 
even. It is known, that the m = states of an even number of spins are eigenstates of S 
with an eigenvalue analogous to (ft.24 ) 

Z xxx = (-1) N / 2 - 1 . (3.27) 



This is a consequence of the fact, that any such state can be given as linear combination of 
states built up as products of I triplet and N/2 — I singlet pairs pSfl . The ground state (the 
vacuum) is singlet, thus all the singlet excitations have the same symmetry as the vacuum 
and the triplets have the opposite. Thus we may consider the degenerated states as forming 
normal SU(2) multiplets of the N spin of a the chain, but equally well we may think of 
these states as the above described multiplets of the dressed particles, and all this is a direct 
consequence of the SU(2) symmetry of the XXX Hamiltonian. In the case of our limiting 
theory there is no Hamiltonian, and the symmetry is recognized in an indirect way: the 
SL is taken through states not showing the SU(2) structure, but this symmetry appears, 
as degenerations characteristic for this symmetry develop. Now we see, that the reflection 
symmetry properties complete this structure. 



G. Boundary conditions in the limiting theory 



In the previous part of the Section we supposed, that the N — > oo limit is taken through 
integer multiples of four (N = 4/V). Actually in order to be able to define the limit of the 
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momentum uniquely and to obtain the parameters in the limiting equations together with 
the properties of the excited states properly defined, we have to fix the way of the N — > oo 
limit, but we have four different possibilities: N = AM + u, M — > oo and v = 0, 1, 2, 3. In 
these cases ( |2.15|) and (|3.24j ) are 



Q{I , N) = Q(I Q , u) = 27r(/ - 1/2 + u/A) , (3.28) 

E(/ 0) N, I) = E(J , u, I) = (_i)(2/o-D + ^ (3 . 29) 

(the latter applying for v — 0, 2 only, as £ has no eigenstates among the BA states of the 
iV=odd chain). Now we want to investigate the v = 1, 2, 3 cases. 

First consider the z/ = 2 case. Now we have to define the momentum in the continuum limit 

as 

P = g(/ °' 2) ~ 7r + £ M sinh(^) (3.30) 
a h 



or 



p= Q(/ ,2) + ^ MsinhW ( 331 ) 



to pick the J = 1/2 or I = states respectively. In this procedure all the results are exactly 
the same as before, except those concerning the reflection symmetry (and as we shall see the 



S- matrices must be affected too). Now, according to (|3.29|) (( [2.16|) ) the secular equations 



( p. 17] ) with I = describe the excitations of the reflection symmetric vacuum, and Iq = 1/2 



gives the excitations of the antisymmetric one (just the opposite as for v = 0). The nature 
of the excitations should be given by the symmetry of the corresponding vacuum, and the 
secular equations are the quantizations of the momenta of the individual particles. Now 
the fact, that the same set of particles (say a singlet pair above one of the vacua) may 
obey two different quantization condition, can be interpreted so that in the two cases the 
boundary conditions are different: if they are periodic in one case, they are anti-periodic in 
the other (the difference in the twisting is it). The two particle scattering matrices given by 
( |3.23j ) are obtained supposing periodic boundary conditions. From the very same secular 
equations but supposing anti-periodic boundary conditions the same S'-matrices multiplied 
by (—1) are obtained, i.e. interpreting the limiting theories for v — 0, Jo = 0(1/2) and v = 2, 
Io = 1/2(0) as being the same ones but with different boundary conditions, is consistent with 
the S'-matrices obtained. It is also consistent with the reflection symmetry of the 5* z = 
states: we checked numerically (for iV = 18), that also in the v = 2 case the two particle 
states scaling to singlets have the same symmetry as the vacuum. 

Now it is not hard to see, that with slight modifications of the prescription of the limiting 
process for the momentum we can consistently define the scaling limit also for v = 1 and 
v = 3. This way, together with the two possible values of Jq we can obtain four cases. 
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Also these are described by the equations ( |3.17H3.19|) but with the restriction, that n(i?) 



must be odd. Also in these cases the eigenstates form SU(2) multiplets corresponding to 
the different combinations of n($) spins of length 1/2 and Eqs. ( |3. 17^3.1^ ), just as in the 
previously discussed cases, give directly the I = m = n($)/2 — n(n) states. We interpret 
these four sets of states as the excited states with odd numbers of particles above the two 
vacua, both with two different boundary conditions. Unfortunately, as none of these sets 
contain the corresponding vacuum, we can not tell, which values of the parameters Iq and v 
correspond to which vacuum. Now, however, the boundary conditions can be seen directly: 
the equations of the one particle states read 

Lp($) = 2tt(7o + 1/4 + integer) (3.32) 



indicating, that these particles obey twisted boundary conditions with twisting angle 27t(Iq + 



1/4) |20|J2T 



A possible explanation of the fact that for different z/s we get different limits can be the 
following. In many cases in order to define the local operators in a naive continuum limit 
one has to group the lattice points into elementary cells of certain length P3|JTI1|JTT[ ] , which 
will be the points of the continuum. If the chain length is not an integer multiple of the 
size of the elementary cell, necessarily boundary terms occur in the limit which can cause 
a rapidity independent shift in the phase (twisting) as a particle is taken around the ring. 
Now it seems, that the size of the elementary cell is four, and so there are four different 
possibilities as far as the 'surface terms' are concerned. In addition to this the number 
of particles must have the same parity as the chain-length (the vacua are spinles, while 
the particles are spin 1/2 particles). These together add up to give the variety of sectors 
discussed above. 
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APPENDIX A: 



In this appendix we list those definitions and properties of the elliptic integrals and functions 
| H| we use in the bulk of the paper. 



fcl 



o Vl - k 2 sin 2 a 

K = 2 (1 - k 2 sin 2 a)~^ da 
Jo 

TT 

K' = [ 2 (1 -k' 2 sin 2 a)~^ da 
Jo 



k 2 + k' 2 = 1 

(1 - k') 2 



k 2 



K< = ^(l + \k' 2 + 0(k'*)) 

K = logi + i(logl-l)A;' 2 + 0(A;' 4 - 

K' 7T ^ k 
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4(l + 7T^X + ^(^°)) 



K 21og± v 41og£ 
IT 

7 = 7T 



,2 



k! = 4e 2 ^ — 16e 2 ^ + . . . 

TT 2 2tt 2 _^ 
if = — + e 7 + ... 

2 7 7 

X' = ^(l + 4e~^ + ...) 



dn(u, k) 

cn(w, /c)dn(w, fc) 

— sn(w, fc)dn(u, k) 



d &m(u, k) 

du 
dsn(u, k) 

du 
dcn(u, k) 

du 

ddn(u,k) , 2 , , . 

= — k sniu, k)cniu, k) 

du 

sn(u,k) = sin &m(u, k) 

dn 2 (u,k) = l — k 2 sn 2 (u,k) 

dn(±K + u,k) = k'- ' 



cn(±K + u,k) = Tk 



dn(u, k) 
, sn(u, k) 



dn(u, k) 



F(«^)^logtan(J + ^ 
sn{u) > tanh(w) 
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APPENDIX B: 



Here we analyze the behavior of the vacuum rapidity sets in the SL and give a derivation of 
the energy difference between the two vacua. 

First one can see, that if 7 — > at fixed N the number of rapidities is fixed and they condense 
at the ends of the interval {0, 7r}, (except the A = tt/2 element in one of the sets), and the 
interior of the interval becomes empty. In this limit if the rapidities are rescaled x = A/7 
(if A < 7r/2) and x = (A — tt)/j (if A > tt/2), the rapidity sets of the singlet ground-state 
resp. the lowest energy triplet state of the XXX chain are recovered. If however, the 7^0 



limit is performed under the condition (1.4), although most of the rapidities condense at 
and 7r but the interior of the interval will not become empty, and the properly scaled 
rapidities have a finite density at tt/2, as discussed in [TirB|. This is an important condition 



for the applicability of the method deviced by @ to calculate the energy difference of the 
two vacua. 

The expression for the energy difference of the two vacua given in [|7]] gives the leading terms 
of a double series of the powers of 1/y/N and k^ 2 in which the coefficients are functions 
of 7. If 7 — >• at fixed iV this expansion breaks down as fci — > 1 and also the expansion 
coefficients may explode. (While the expansion breaks down, the true energy difference will 
be oc 1/iV.) For this reason one has to be very careful, and if our (|1.3| - |1.4| ) limit is performed, 
the limit is to be calculated in all terms. Evaluating the terms presented in |7j indicates that 
the limiting energy difference is indeed a function of L (just as (1.5) is). In the following 



we present a simple derivation of ( |1.5| ) in which one can see, that the neglected terms decay 
faster. 

According to the energy difference is given by 

n+ij/2 



AE ~ -4Re / e(9)a (6)e lNq{e) d9 (Bl) 



0+17/2 



where e(9), cr (#) and q{6) are given by ( |2.21 ), (2.12) and Q2.22| ) respectively, and the ne- 
glected terms are of the type 



I 

0+17/2 



j f(6)e iNnq{e) dd (n>2 integer) . (B2) 



The function q{6 -M7/2) has an extremum at 9 = tt/2, and we can use the saddle-point 
method. Introducing the variable $ of (|3.3| ) we find that around 6 = tt/2 

q{6 + ^ 1 /2) = \ In (±^\ + ^ . (B3) 
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i.e. (as k'/a -> M): 

iNq(6 + 11 /2) = N In P^") ~ ■ ( B4 ) 

It is clear, that for large L to pick up the leading contribution we may use this expression 
in the integral, and we may expand also the e and a in terms of This leads to 

A^ 8A ™ ,)2 fl^gf /W^M (B5) 

(27T \ K I J 

v 7 



which yields ( |1.5| ) indeed. We can also see, that the neglected terms ( |B2| ) are 

oc e~ nLM (B6) 

i.e. for large enough size L the energy difference decays as ( |1.5|) . 



APPENDIX C: 



In the bulk of the paper we derived the secular equations of the limiting model as the 
limit of the secular equations of the lattice model. This procedure is meaningful only if the 
approximations used to derive the original equations remain valid in the SL too. 

One of the approximations applied is the replacement of summations by integrals. The error 
introduced this way can be calculated by the method of 0j, and as it leads to the evaluation 
of formulae of the type (|B~I|), we may conclude, that this type of corrections are exponentially 
small in LM. 



Another subtle point of the derivation of the higher level BA equations is the representation 
of each close pair by a single number, namely by its center. As it is known, the close pairs of 
the original rapidity set are defined by those x solutions of the higher level BA equations, for 



which Im% < 7/2 through the rule ( 2.20a ). This approximation is allowable, if exp{— iVfi} 
will remain small in the SL. Now we show, that this is indeed the case. Actually according 

to § 



In 



Z7 



-MX*+ 7) 



(CI) 



with 



E 



?!/) 



{cosh - (x + rmr) + sin 
1 7 7 
cosh ^(x + mTC) — sin(^y) 



(z+ = x + iy, y<i). (C2) 



In studying the behavior of the Q one has to distinguish three cases (see Appendix 0): 
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1. The x scales to or ir: 

, Pn (03) 

In this also y oc 7, f2 —> finite, and NQ — > 00, i.e. z + — z — 17 — > exponentially 

fast in iV. 

2. Another possibility is, that x is finite but a; -f-> ir/2. It is not hard to see, that in this 



case 



jVn ~ ^£ef>5-*> mn(-y) . (C4) 

Also this expression diverges as 7 — > 0. 

3. Finally, if x scales to tt/2, i.e. x = 7r/2 + 7ReK/7r (and y = jlnm/ir + 7/2 with 
ImK < 7r/2), one finds, that 

NQ = LM cosh(ReK) cos(ImK) . (C5) 



The close pairs described as case 1. and 2. form exact 2- strings or quartets in the scaling 
limit, but they disappear from the secular equations. Nevertheless they play an important 
role in determining the multiplet structure of the states. The close pairs of case 3. are in 
strong analogy with the close pairs of a finite lattice model in the sense, that they form 
exact 2-strings or quartets in the infinite size limit only. As, however, the deviation is 
exponentially small in the size of the system, the representation of the close pairs by their 
centers is justified. 



APPENDIX D: 



In this Appendix we discus the structure of the solutions of fl2.18| ), what we write in the 
form 



ffi sin( Xi - 9 h + ry/2) 
W sm( Xj -9 h - i 7 /2) 



n jT sin(xj - Xk + in/) 
^ sin(xj -Xk~ *7) 



(Dl) 



We give an account of their different types, and relate their 7^0 limit to the solutions of 



n 



Kj — $h + *7r/2 



\ Kj -"&h- in/2 



7l(ft) 

n 



k=i 



(D2) 
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In this analysis we shall keep the scaled rapidities fixed while seeking for solutions to Xj 
as a perturbative expression in 7: 

0& = £ + -0a, Xi = X,0 + -^ + O(7 2 )- (D3) 

2 7T 7T 



Note that Xjo — f and Xj = Kj for the rapidities which scale to tt/2 as was defined in ( |3.12 ). 



We label the types of solutions by three numbers {H, I, m}. In this label the first number is 
H = n(i?), the number of particles, the last one is m = S z = n($)/2 — n(x), i.e. the spin of 
the corresponding state, and the middle one is given by I = n($)/2 — n(n). This quantum 
number is specified in the bulk of the paper, when the multiplet structure is discussed. 
(Note, that I — m = n(x) — ti(k) > is the number of x s n °t scaling to tt/2.) 

We analyze in detail the states with n(#) being equal to two and four. 

States with n(?9) = 2 
n(x) = 

The simplest such state is the one, in which there is no x- The label of this state is {2, 1, 1}. 
n(x) = 1 



The equation 



has two solutions: 



sin( X - 9 h + » 7 /2) = 



X = — 2 — ' 

X = — — (mod7rj . (Do) 



In the first x scales to 71/ 2: 



with 



X = ^ + \ (D7) 

7T 2 



/t = , (D8) 
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in the second not. According to our notations the two solutions are labeled by {2, 0, 0} and 
{2, 1, 0}, respectively. 

Note, that in the SL the {2, 1,1} and {2, 1, 0} solutions pairwise become degenerate, as the 
equations fixing the tf^'s will be exactly the same for them. 

States with n($) = 4 

n{x) = 

The simplest state is the one with no x- This is labeled by {4, 2, 2}. 

n{x) = 1 

The states with one x are characterized by the four nonequivalent solutions of the equation 

L\ MX -0 h - < 7 /2) • lUyJ 

Three of these are well approximated (in the 7^0 limit are exactly given) by the three 
(real and finite) solutions of the equation 

and one x is given in leading order by 

x=\h*-\ ( m ° d7r ) = 1 e *h ( m ° d7r ) • ( Dn ) 

h=l h=l 

The first three solutions have the label {4, 1, 1}, while the fourth is {4, 2, 1}. 

nix) = 2 

The solutions for the equations of two x s 

^ sin(xi -0 h - i 7 /2) ^ sin(x; - Xk - n) 

are of three type. The set of first type coincides in the 7 — > limit with the set of solutions 
of the equations 
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M (D13) 



i.e. they are given through (P7|) by the solutions of 



Xj 




— 


of 






K i 


- Kfc 


+ «7T 


« 3 - 


- «fc 


— in 



■q K j ^fe + vkJ2_ _ j-j ,^ ,v fc i ». (D14) 



These solutions are labeled by {4,0,0}. (These solutions can be given even analytically. 
The two nonsingular solutions (i.e. in which the roots are not at the singularities of the 
equations) are of the form 

4 ft 

with (A$) 2 given by 

((Atf) 2 ^ 2 = 1 (-23 ± V# 2 + 12C) (D16) 



with 



£ = 2 U + EW(<^0, c=(-) -2 - (Di7) 

vz/ ft<ft' vz/ vz/ ft<ft' ft 

and 

^ = ^-jEW- (D18) 
4 ft' 

(The singular solutions are non-physical ones as they lead to vanishing wave-function.)) 
Another set of solutions is given by 

7ft 71 

Xl = V + 9 ' 
7T Z 



X2 2 



i ^ - 2 Xi) - \ (modvr)-^ #ft - 2ft) (modTr) , (D19) 



with ft being one of the solutions of ( piO| ). The label of these three solutions is {4, 1,0}. 
Finally there is one solution labeled by {4, 2, 0}, for which 

Xi,2 = \ E h ± i tanh" 1 (l/y/3) + 0(j) . (D20) 
4 ft 
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The {4, 2, 2}, {4, 2, 1} and {4, 2, 0} solutions become degenerate in the SL as the equations 
determining their fi^s become identical. For the same reason the {4,1,1} and {4,1,0} 
solutions form degenerate pairs having the same k. 

The general case 

Now we attempt to classify the solutions of the equations ( pip in the (n(9) > 2n(x)) general 
case. It is clear, that in the scaling limit in each solution some of the x s (say < n(x) 
of them) will scale to tc/2 in the form ( D7 ) with the ks satisfying ( |D2| ) (with n($) = n(6)). 



Now our task is to give the number of solutions at given n(9), n(x) and tz(k), i.e. the number 
of solutions labeled by the same label {H, I, m}. Giving this number we have to make two 
points. First we suppose, that an equation (DTI) has 



Nxxz(n(9),n{x))= (D21) 

nonequivalent solutions (one solution being a set of xs) . This statement is strongly related to 
the question if the BA yields a complete set of eigenstates for the easy axis XXZ Heisenberg 
chain. Our other assumption is, that an equation of the type ( p2|) has 





N X xx(n(#),n(K)) = T - T\ (D22) 



different solutions. This is related to the question of completeness of the BA solutions for the 
isotropic case (in which a subset of the HLBAE are again of the type (H2). The difference 



of the two statements originates from the fact, that in the isotropic case the BA gives the 
S 2 = S Z (S Z + 1) states only (i.e. not all of the states), while in the anisotropic case S 2 is not 
a good quantum number, and the BA should give all the S z > states. Now our conjecture 
is, that the number of those solution of ( pi| ) in which ti(k)(< n(x)) of the x s scale to 7r/2 
in the scaling limit is again (|D22| ) (with n{9) = n($)), and in each solution the n(n) of the 



Xs which scale to n/2 are given by the different sets of ks solving ( p2[ ). This means, that 
if n(x)(> n(K,)) and a k set solving ( P2|) are given, there is no further freedom, and there 
is only one possibility to complete the set of xs generated by the ks to a set solving ( pl[) . 
In other words, in a class of solutions {H, I, m} there are Nxxx{H, \H — I) different states 
uniquely determined by the different solutions of (|D2| ) independently of the value of m, (for 
which the only requirement is I > m > 0). This can be seen in the cases discussed, and is 
supported by the fact, that 

E N X xx(n(e) : n(/c)) = N XX z(n(6), n( X )) ■ (D23) 

n(re)=0 



A consequence of this is that in the scaling limit the solutions can be collected into groups 
of degenerate states: at a given H and I the members of such a group are characterized by 
common dh and sets but differ in m what runs from to I. 
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APPENDIX E: 



In this appendix we show some solutions of the BAE of the finite lattice which has S z = 
and we give the symmetries of these states under E. In these examples the length of the 
chain is always 20, the number of holes is 2, and we have one close pair. The equation 
(2.8) was solved by iteration which was started from the vicinity of the solutions of the 
HLBAE ( |2.17| - |2~T8|) . After this the amplitudes ( |2.2j ) of some configurations (which we think 
dominate) were calculated and it was found that any of these amplitudes was ±1 times the 
amplitude of the corresponding reflected configuration. (The relative error was 10" 6 but it 
is plausible that the reflection symmetry holds exactly and the error comes only from the 
numerics.) The iteration converges fast also for much larger N, but the time needed for the 
calculation of the amplitudes grows rapidly with N resulting a strong upper limit for the 
number of lattice sites we can handle. For illustration, Table | shows some of the solutions 
we have found. There are two types of solutions: in type 1 the real part of the close pair v c 
is between the two holes, while in type 2 it is located outside. It can be seen in Table | that 
in case of Iq = 1/2 type 1 states are symmetric and type 2 ones are antisymmetric under E 
while the opposite is true for I — 0. It should be noted that the connection between the 
type, Iq and E is independent of the interval the rapidities are collected in because shifting 
one of the 9s by it both the position of x compared to the holes and the value of Iq changes 
|| . In our choice (0 < v < tt) type 1 resp. type 2 solutions scale to singlet ( |D5|) resp. triplet 
( |D6| ) states as it is indicated in Table [I]. 
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TABLES 



T 
1 


(7, U 


7, v c , I , E 




r 
i 


(/, u 


T V 


-0.5 


0.0462991 








-1.5 


0.1438185 




-1.5 


0.1481864 




7 = 1 




-2.5 


0.3092072 


7 = 1.5 


-2.5 


0.2663213 








-3.5 


0.5164954 




-o.o 


(J. 4241941 


V c 


= 1.5046086 ± i0.6433822 




-4.5 


0.8408343 


v c = 1.0670662 ± i0. 7483516 


-4.5 


r~i r~7f\ A C\ A A'l 

0.7049447 




(— » singlet) 




-5.5 


1.6284527 


(— > singlet J 


-5.5 


2.4271298 








-6.5 


2.2989503 




-6.5 


2.7097631 




J = 2 




-7.5 


2.6116381 


I Q = 


-7.5 


2.8692656 








-8.5 


2.8192199 




-8.5 


2.9880692 




E = +1 




-9.5 


2.9854038 


V 1 1 

2j = —1 


-9.5 


3.0904861 








-10.5 


3.1355551 




-0.5 


0.1361835 








0.5 


0.0407303 




-1.5 


0.3022448 




7 = 1.5 




-0.5 


0.1421697 


7 = 1 


-2.5 


0.5124295 








-1.5 


0.2594036 




-3.5 


0.8528326 


Vc 


= 0.8744719 ± i0.7496363 




-2.5 


0.4157572 


v c = 1.3736073 ± z0.6398024 


-4.5 


1.6684991 




(— > triplet) 




-3.5 


0.6965562 


(— » triplet) 


-5.5 


2.2976980 








-4.5 


2.4190131 




-6.5 


2.6042345 




4 = 2 




-5.5 


2.7023676 


io = 


-7.5 


2.8110406 








-6.5 


2.8629747 




-8.5 


2.9773204 




E = -1 




-7.5 


2.9825513 


E = +1 


-9.5 


3.1276728 








-8.5 


3.0850192 





TABLE I. Some solutions of the BAE flU) of the N = 20 lattice with 8 real rapidities, 2 
holes and one closed pair. In the first two columns the real rapidities v, the holes 9 (emphasized 
by bold face) and the corresponding quantum numbers / can be found. The third column shows 
the parameter of the anisotropy 7, the rapidities of the closed pair v c and the value of Iq and E. 
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